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A new method for investigating the dynamics of atomic magnetic moments in current-carrying 
magnetic point contacts under bias is presented. This combines the non-equilibrium Green's function 
(NEGF) method for evaluating the current and the charge density with a description of the dynam- 
ics of the magnetization in terms of quasistatic thermally-activated transitions between stationary 
configurations. This method is then implemented in a tight-binding (TB) model with parameters 
chosen to simulate the main features of the electronic structures of magnetic transition metals. We 
investigate the domain wall (DW) migration in magnetic monoatomic chains sandwiched between 
magnetic leads, and for realistic parameters find that collinear arrangement of the magnetic mo- 
ments of the chain is always favorable. Several stationary magnetic configurations are identified, 
corresponding to a different number of Bloch walls in the chain and to a different current. The 
relative stability of these configurations depends on the geometrical details of the junction and on 
the bias, however we predict transitions between different configurations with activation barriers of 
the order of a few tens of meV. Since different magnetic configurations are associated to different 
resistances, this suggests an intrinsic random telegraph noise at microwave frequencies in the I-V 
curves of magnetic atomic point contacts at room temperature. Finally, we investigate whether or 
not current induced torques are conservative. 

PACS numbers: 75.75.+a, 73.63.Rt, 75.60.Jk, 72.70.+m 



I. INTRODUCTION 

Most of the intriguing properties of ferromagnetic 
nanoscale atomic structures arise from the close inter- 
play between magnetic phenomena and electronic trans- 
port. As the magnetization can be controlled on a length 
scale smaller that the spin diffusion length of the conduc- 
tion electrons^, the spin scattering is affecting the over- 
all resistance of an atomic ferromagnetic device. This 
is the principle behind the giant magnetoresistance ef- 
fect (GMR)2i£. Remarkably also the opposite effect is 
possible, i.e. the electronic current, as proposed by 
Slonczewski^, can transfer spin and alter the magnetic 
configuration of the underlying ferromagnetic structure. 
Magnetization switching, caused by spin-polarized cur- 
rents, has been observed experimentally in point contact 
measurements^, and in nanopillars 7 . 

It is then clear that the modeling of these atomic-scale 
ferromagnetic devices requires the combined description 
of electronic transport and of the magnetization dynam- 
ics at the atomic level. For this purpose we have devel- 
oped a general scheme for evaluating spin-currents and 
associated current-induced torques, which allows us to in- 
vestigate the magnetization dynamics and the transport 
of magnetic point contacts under bias. Our problem and 
our method mimic closely, in philosophy, electromigra- 
tion problems (thermally activated current-driven struc- 
tural rearrangements), where now the direction of the 
local magnetic moments takes the place of the atomic 
positions as "reaction coordinate" . 

Although our scheme is general and is conceptu- 
ally transferable to first-principles Hamiltonians (for in- 



stance, within density functional theory), here we apply 
the method to a simple self-consistent tight-binding (TB) 
model. This has the benefit of being reasonably realistic 
while keeping the computational overheads to a mini- 
mum. 

The paper is organized as follows: in the next section 
we introduce our method, describe the model and sketch 
out the techniques, used for the calculations. Then, in 
section II. C, we discuss our approach to the interplay 
between spin-polarized transport and the magnetic con- 
figuration, applied to a particular atomic structure. In 
section III we report a set of results, which explore the 
stability and the activation barriers for transitions be- 
tween various MM configurations under bias, as well as 
the effect of the model parameters on the physical proper- 
ties of the system. Finally, we carry out a numerical test 
to see whether or not the torques in these open-boundary 
non-equilibrium systems are conservative. 



II. THE METHOD 

A. General Idea 

Our scheme for studying current-induced dynamical ef- 
fects of the magnetization in atomic-sized nanostructures 
is a generalization of the combined quantum-classical 
dynamical methods used in electromigration problems 8 . 
Here we treat the magnetic degrees of freedom as clas- 
sical variables and the conduction electrons as a quan- 
tum system. This is appropriate when the magnetic mo- 
ment (MM) arises from some deep orbital levels, such as 
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in the case of rare earth ferromagnets, but it may ap- 
pear questionable for magnetic transition metals (Fe, Co 
and Ni), where the d electrons responsible for the mo- 
ment also take part in the conduction™. However, since 
the Coulomb energy is orders of magnitude larger than 
any energies connected with the electron flow, it is safe 
to assume that only the direction of the local atomic 
MM is affected by the current but not its magnitude. 
This effectively is an adiabatic approximation, in the 
spirit of the Born-Oppenheimer approximation for the 
nuclear dynamics, where now the orientation of the lo- 
cal MMs is a slow variable compared with the internal 
electron-electron interactions that set the magnitude of 
the MMsifi*ii. The Hamiltonian for the combined con- 
duction electron/MM system can be then written in gen- 
eral as 



H({<l>})=H e + V es ({<P}), 



(1) 



where we have isolated the interaction of the conduction 
electrons with the local MM, V cs ({cf>}), from the "free" 
electron Hamiltonian H c . In this framework the local 
moments are uniquely specified by a set of angles {</>} 
with respect to a given direction. 

We may now write down the generalized forces (in this 
case, torques) conjugate to the classical variables {<p}: 



(2) 



where 1^) is a state vector of the electronic system. 
Equation (0) has the appearance of the usual Hcllmann- 
Feynman theorem for stationary states. However it is 
valid in general dynamical situations, for systems driven 
arbitrarily far from equilibriumi2ii£. 

The set of equations Q and J5J, combined with an 
appropriate method for calculating the non-equilibrium 
electron state vector |\&), and therefore the current, is 
the basis for our method for describing the interplay be- 
tween transport and magnetic properties. In this work, 
we seek to map out the activation energy barriers for 
magnetic rearrangements, in order to determine the pref- 
erential magnetic configurations of the system and to 
study transitions between them. We achieve that as 
follows. First, we seek the stable configurations. We 
evaluate the non-equilibrium state vector \^>), in a one- 
electron picture, for a given starting MM configuration 
{(f) } by solving the scattering problem associated with 
the Hamiltonian H({<fi a }). Then, by using equation @ 
the torques for that configuration are calculated. Then 
static relaxation of the torques is carried out (by steepest 
descent), while recalculating the self-consistent current- 
carrying electronic structure, and torques, at every step. 
The procedure continues until all torques vanish. This 
condition gives the stationary magnetic configurations 
{<&}. 

Once the stationary magnetic configurations have been 
found, we can calculate the activation energy barriers 
for thermally activated transitions between two different 



configurations {c[> lmtla1 } an d {$ fina1 }. We choose a single 
classical dynamical variable (f>j as the reaction coordinate 
and rotate it its starting value <^ mtlal to its final value 
^finai^ every step on the way the torques acting on 
all other MMs are kept relaxed to zero. The work done 
by the classical degrees of freedom during this quasistatic 
transition is then obtained by integrating the torque on 
the reaction coordinate <pj over the migration path. The 
work done over the full transition is 



W 



TV 



(3) 



The energy barrier profile, on the other hand, is given by 



(4) 



where = {<j>x{4>j), foifa), ••, fa, — , $N(fa)} is the 
magnetic configuration, for a given <j)j, defined by the 
condition Ti = for every i =/= j. 



B. Transport Method 

The Keldysh non-equilibrium Green's function 
(NEGF) method is used here for describing the 
transporti^ii^. We expand the Hamiltonian H({<j>}) and 
the one-electron wavefunctions in a localized atomic 
orbital basis set, and we decompose our system into two 
current/ voltage probes sandwiching a central region. 
For an open current-carrying system H can be written 
as 



H = Hr, + Hn + Hr 



H 



LC 



H 



RC 



(5) 



where we have introduced the Hamiltonian for the left- 
(right-) hand side current /voltage probe H^ = H^({(j)}) 
(Hji = ,{{4>})), that for a central scattering region 
He = Hc({0}), and the coupling matrix between the 
left (right) contact and the scattering region H^c = 
H L c(M) (H RC = iJ RC ({<£})). The latter are indeed 
spin-polarized operators, i.e. H = X)o-H CT , but for the 
sake of simplicity of the expression we skip the index a 
in the following formulas, thus refering to either of the 
spin components. 

The NEGF method allows us to map this in principle 
infinite problem (Hr and Hl are infinite matrices) on 
an auxiliary finite problem. The key observation is that 
one can describe the effects of the current / voltage probes 
over the scattering region by means of their correspond- 
ing self-energies Sl and £r. These are non-hermitian 
matrices which contain all information about the elec- 
tronic structure of the probes and their occupation. They 
can be written as 

s l = Hl c g L H LC , and S R = H r c5rH rc , (6) 

where we have defined the surface Green's function (GF) 
for the left- (right-) hand side probe <7r (<7l)- Hence the 
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"effective Hamiltonian" of the scattering region in the 
presence of the current/ voltage electrodes is written as 

ffcff = H c + Sl + Sr. • (7) 

Note that this is a finite non-hermitian matrix. Conse- 
quently the number of electrons in the scattering region 
is not conserved. 

Now we can construct the retarded GF associated with 
the scattering region plus the leads 

G(E) = lim [(E + »C) - H c - S L - Er]" 1 , (8) 

and the associated (non-equilibrium) density matrix 
1 f°° 

P=- dE (n L (E)r] F (E) + n R (E) V K(E)) , (9) 

^ J — OC 

where n A (E) = G(E)T A (E)G jf {E) are the partial den- 
sity of state operators for electrons originating from each 
lead (A = R,L), T = i(E(E) - Y^(E))/2 is the non- 
Hermitian part of the self-energy, T]p{E) = T)p(E — p A ) 
are the corresponding Fermi distribution function for the 
electron reservoirs with chemical potential \x a- 

Under our basic assumption of "reflectionless" leads 
we can decouple the subsystems of electrons originating 
from the left and the right lead and treat them as sep- 
arate statistical systems with electrochemical potentials 
Pl(r)- In equilibrium /i^ = pn = p, and a finite bias V 
is introduced as Pl(r) = £t± \e\V/2, so that the electron 
flow is from the left to the right lead for V > 0. Thus the 
bias V is assumed not to change the electronic structure 
of the leads, but only to rescale the energy levels. In prac- 
tice V is introduced as a rigid shift of lead-Hamiltonian 
on-site energies 

Hl/r — H h/R ± -I, (10) 

where X is the identity matrix for the respective lead. 

Self-consistency in our calculation is introduced by as- 
suming that for a given magnetic configuration {0} the 
Hamiltonian He depends solely on the scattering region 
density matrix He — Hc[p]- This is equivalent to as- 
suming that the underlying electronic structure theory is 
a density-based theory, such as Hartree-Fock or density 
functional theory. In this case the set of equations J7J), 
© and 10 defines the self-consistent procedure. First 
one computes the scattering region GF (equation (JSJ) for 
Hc\fio] evaluated at an initial density matrix pg. Then 
from the GF a new charge density p\ is calculated and 
used to construct the new Hamiltonian He [pi] ■ This pro- 
cedure is iterated until reaching self-consistency, that is, 
until p n+ i = p n . 

Finally from the converged GF the net current is cal- 
culated asi2J£ 

p r°° 

1{ y) = jll / dE (Vf(E)-Vf'(E)) 

Tv(T L (E)G^T R (E)G(E)) (a) (11) 
where we have summed over the spin index a. 



C. The Model 

The techniques described in the previous sections are 
general and can be applied to a large class of Hamiltoni- 
ans. In this work we focus our attention on a simplified 
model, which contains the fundamental ingredients for 
describing a current-carrying magnetic point contact, but 
at the same time does not present massive computational 
overheads. The structure we investigate is schematically 
represented in figure ^ It consists of two semi-infinite 
leads with a simple cubic lattice structure and a 3 x 3- 
atom cross section connected through a linear chain of 
three atoms. Each atom carries a local magnetic mo- 
ment, our classical quantities, arising from the deeply 
localized (i-electrons. The magnetic configuration of the 
leads is fully polarized (all MMs in a given lead point 
in the same direction) and we investigate the situation 
where the magnetizations of the two leads are opposite 
to each other. In contrast the three MMs of the atoms in 
the chain are allowed to rotate. A given magnetic config- 
uration of the chain is thus described by the three angular 
coordinates (<^>i, 4>2, ^3), with the convention that <f> = 
{<f> = 7r) for a spin alignment parallel to that of the mag- 
netic moment of the left (right) lead. We consider spin 
rotations only in the x-y plane thus neglecting the longi- 
tudinal angle (as in a Bloch wall). The alternative choice 
is to consider MM rotations in the x — z plane (a Neel 
wall), but as far as we neglect the magnetocrystalline 
anisotropy, these two models are identical. 




FIG. 1: Scheme of the magnetic point contact 



The current is carried by electrons belonging to an 
s band, which is described by means of a single-orbital 
(plus spin) TB model. The Hamiltonian of equation JTJ 
is therefore explicitly written as 

H(W) = X)[0ffe)ij + (^int)y] c\ Cj + Kla SS ical({</>}) , 

(12) 

where c\ and Cj are creation and annihilation operators 
for electrons at the atomic sites i and j respectively. The 
matrix elements of the "free" electron part are those of 
a nearest-neighbor TB model 

(ff e )« = [co + U {pi - P °)]5 l3 + 1 5 l J±1 , (13) 
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where eo is the on-site energy, 7 is the hopping parameter, 
U is the on-site Coulomb repulsion, p° the reference on- 
site charge corresponding to the neutral free atom, and 
p is the self-consistent local charge. The potential V es 
of equation Q has now been separated in two parts: 
Vi n t and ^classical- The interaction between conduction 
electrons and the local MMs is contained in Vi n t which in 
our model reads 

(Vint) ij = ~ s i~2 cos< ^ ; 5i i ■ ( 14 ) 

T I 

where Sj = p\ — p\ is the local spin polarization of the 
s electrons at site i and J is the exchange parameter. 
Therefore V- m t describes a Heisenberg-type interaction 
between the local classical MMs and the current carry- 
ing s electrons. Finally, the classical term T4i ass icai pa- 
rameterizes the interaction between local MMs. Here we 
assume a Heisenberg spin-spin interaction 

^classical = ~ ~7T & ' &j = ~ X/ C ° S ^ i ~ ^ ' 

(15) 

where Jdd is the intersite exchange integral and we have 
assumed normalized classical spin \S%\ = 1, in such a way 
that \8i\ is incorporated in the definitions of J and Jdd- 
In summary our model is that of s conduction electrons 
exchanged coupled to local magnetic moments, in turn 
described by a Heisenberg-type energy. This is usually 
known as the s-d model^. 

The torque experienced by the i-th. local MM in the 
chain is then obtained from equation (J2J and reads 

Ti = -^-Si sin (j)i - -^(sin(0j - + sin(0 i+1 - fc)) 

(16) 

where i = 1,2,3, and we have defined 0o = 0, 04 = tt 
since the magnetization of the two leads is considered 
pinned in an antiparallel arrangement. 

In this simple model the surface GF (at a general com- 
plex energy E) of the leads have an analytical form. In 
reciprocal space, 

9 (E,k) = E -^-V[-p ^^, Im(E)>0(17) 

where e(k) is the energy, as a function of transverse 
wavevector (in appropriate unites) k = (k x ,k y ) with 
k x = 1 , . . . , N x , k y = 1 , . . . , N y for an N x x JV^-atom simple- 
square monatomic slab in a nearest-neighbor orthogonal 
TB s-band model, 

c(t)=eo _ 2 | 7 | COS ( ] ii_)_ 2 | 7 | C08 (_^_). (18) 

The expression of equation 117fl is then expanded over the 
real-space basis^i and used in the matrix equation for the 
self-energies. The definition of the complex square-root 
is given iniS. 



III. RESULTS 

Here we investigate the magneto-dynamics of atomic 
point contacts, and in particular of the model structure 
described in figure^] The TB parameters are eo = —3 eV, 
7 = — leV, U = 12 eV, which gives a large bandwidth 
for the s electrons and provides local charge neutrality 
as expected in a metal. For the exchange parameters 
we investigate the range < J < 3 eV and < Jdd < 
5 eV. However we have identified the values J = 1 eV and 
Jdd = 50 meV as a realistic choice for simulating the main 
physics of magnetic transition metalsiSii2iS£, and we will 
refer to those values as the "realistic parameters" . 

We start our analysis with studying DW migration in 
the three atom chain. As in figure ^ the magnetic mo- 
ments of the leads are in the antiparallel configuration, so 
that that a DW nucleates in the chain. We then investi- 
gate the displacement of the Bloch wall from the interface 
between the first and the second atom in the chain to that 
between the second and the third (generated by a rota- 
tion of the magnetic moment of the middle atom) . These 
simulations use the realistic parameters given above, so 
they can be related to point contact experiments^-. Both 
the cases of spatially symmetric and asymmetric chains 
are studied. Then we explore the effect of varying the 
strength of the exchange parameters, and we identify 
three different regimes. Finally, we revisit the problem of 
whether or not generalized forces away from equilibrium 
are conservative, and demonstrate numerically that the 
torques in the present system under current flow are not 
conservative. 



A. Domain wall migration 

By performing numerical minimization of all the 
torques, exerted on the MM in the constriction, with 
various initial conditions, we have determined that all 
eight collinear arrangements, such as (0,0,0), (0,0, n), 
(0,7r,7r), (tt, 7r, 0) etc. are stable zero-torque magnetic 
configurations and we have studied various transitions 
between them (see fig llOH . In particular, we have investi- 
gated in detail the migration of an abrupt DW within the 
atomic chain, i.e. the transition between the (0, 0, tt) and 
(0, 7T, 7r) magnetic configurations, achieved by rotation of 
Mi as described above. Physical characteristics of this 
process as function of 02 are presented in Fig[21 It is 
observed that during the rotation of M2 its neighboring 
MMs experience small tilts from the collinear alignment 
and after a turning point fall back to their initial state. 
The intersite exchange coupling is not strong enough to 
induce spin flips of the neighboring MM along with the 
one that is rotated and even hypothetical values of Jdd 
up to 0.4 eV do not change this picture (see figEl- This 
observation suggests that the dynamical processes of the 
MMs in the constriction can be decomposed into series 
of single MM rotations. 

The torque T2, computed as a function of 02 (figGt) 



5 




§Jn (rad) (j> 2 /jt (rad) 



FIG. 2: Typical calculation of microscopic properties during DW 
migration in the contact (for J = leV, Jdd = 50meV) as function 
of the "reaction coordinate" cfe: (a) The stable angular variables 
(j>i and (f>3\ (b) The three on-site spin polarizations s; = p\ — pj; 
(c) Torque and work, performed by the MM; (d) Net current at 
V = 0.5 V. The voltage in panels (a)-(c) is zero. 




<j) /it (rad) (j) J% (rad) 

FIG. 3: Effect of the external bias on the microscopic properties 
of the contact as function of (p2- the zero-torque positions of (a) <j>i 
and (b) </>3'i ( c ) the on-site spin polarizations {.Si = p\ — p\}i—i,2,z', 
(d) the work profile; (e) the total net current and (f) its polarization 
«/ = (/ T -/|)/(/ T +/ i ) 



at every point on the way, is interpolated and integrated 
according to © to determine the effective energy barrier 
for the DW migration 

W(fa) =- / 2 T 2 d^. (19) 
Jo 



Because of the specific geometric and time-reversal 
symmetries of the system, the two states (0,0, tt) and 
(0, tt, 7r) are macroscopically identical. Thus the cal- 
culated energy barrier between them is symmetric and 
the total work W(ir) for the quasi-static process is zero 
(fig[2t)- The activation energy for this process in our 
TB parameterization is 54 meV. It is found that the con- 
ductance of our system depends on the alignment of the 
MM and in this case the net current shows a symmet- 
ric bell-shaped dependence on fa (figEJi). For this case 
(V = 0.5 V), the conductance varies from 1.76e 2 /h at 
the collinear states fa = 0, tt to a maximum of 1.86e 2 //i, 
reached at fa = vr/2. 

Further, it is observed that the external bias, driving a 
spin-polarized current, suppresses the response of fa. 3 to 
the motion of fa (figE^,b) but enhances the onsite po- 
larizations (figEt) as well as the energy barrier (figEli). 
At any finite temperature, this phenomenon would man- 
ifest itself as suppression, with increasing bias, of the 
frequency of DW transitions back and forth between the 
two stable magnetic configurations. The net current pro- 
file is slightly sharpened as the bias increases (figOt) an d 
it also becomes more spin-polarized due to the increased 
misalignment of the correspondent spin-polarized bands 
in the two leads (figGU). 



B. Non-uniform contact 

Current-induced relaxation of the atomic positions can 
break the spatial symmetry in point contacts similar to 
ours& and substantially weaken the stability of these sys- 
tems. To investigate the effect of a small asymmetry in 
the contact geometry on the migration barrier for the 
DW, we map the displacement of the middle atom from 
its symmetric position onto a small variation of hopping 
integrals between the middle atom and its neighbors in 
the chain 

712=7(1 + *), 723 =7(1-*)- (20) 

It results in breaking the symmetry about fa — tt/2 of 
the effective energy barrier observed in all our previous 
calculations (fig'EJ. The total work for the (0,0, tt) — > 
(0, tt, tt) transition is negative thus the internal energy 
of the classical MM is increased. The degeneracy of the 
(0, 0, tt) and (0,tt, tt) state is lifted, as the spatial sym- 
metry, associated with a reflection plane at z — 0, is no 
longer present. When the hopping parameters are altered 
by 5 = 5%, as if the middle atom is slightly shifted to 
the left, the (0, tt, tt) configuration becomes energetically 
preferable, alternatively, 5 = —5% favors the (0, 0, tt) 
state, with all the microscopic properties being invariant 
to a change fa — > (tt — fa) (figBJ. 

We have presented the typical microscopic proper- 
ties in figure |SJ As expected the effective ferromag- 
netic coupling between the MMs is further increased by 
the enhanced electronic coupling. The onsite polariza- 
tions of all the atoms shift almost rigidly as the mid- 
dle atom is brought towards one or the other of the 
leads. The net current shows significant asymmetry 
from the regular bell-shaped dependence on fa and the 
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1 O 0.5 

<t>-/7T Crad) 



FIG. 4: Effect of an asymmetry in the hopping integrals in the 3- 
atom chain with S = ±0.05 on the barriers for (0, 0, 7r) — ► (0, 7r, 7r) 
transition at different bias: (a) V = V; (a) V = 1 V(a); V = 2V. 



more stable configuration is always found to be less con- 
ducting (figEJi). The I-V characteristics of the previ- 
ously degenerate (0,0, it) and (0, 7r, 7t) states is split into 
two branches, whose displacement increases with voltage 
(figEf ) and reaches 10% for V = 2 V. Thus we expect DW 
migrations within the constriction, in the case of small 
deviations from a uniform geometry, to be accompanied 
by random-telegraph-noise-like variations in the net cur- 
rent. The interplay between the current-induced relax- 
ation of the magnetic and mechanical degrees of freedom 
is the subject of work in progress. 
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FIG. 5: Effect of an asymmetry in the hopping integrals with 
S = ±5% on the microscopic properties of the contact (as in figlSI 



C. Mapping out the parameter space 

The torques defined in Ijl6|l depend explicitly on the 
exchange parameters J, Jdd and the balance of these 
two coupling mechanisms determines the spin dynam- 
ics in the constriction. Figure EIa,b,c) describes the 
variation of the equilibrium magnetic properties of the 
(0,0, 7r) — > (0,7r,7r) transition as the intersite exchange 
strength Jdd is varied. The effect of greater Jdd is in- 
creased disalignments of the two neighboring MM <j> 13 
following the rotation of fa (&,b). Because of this the 
corresponding spin-polarizations si t s (but not S2) are 



slightly affected (jBJ;) as a result of the stronger coupling. 
The typical bell-shape of the current vs. fa is broadened 
as Jdd increases (figEfc) due to the fact that stronger in- 
tersite exchange coupling tends to make the three MMs 
in the contact more uniformly distributed in angle, which 
makes the contact better-conducting for any fa ■ Also the 
current is somewhat less polarized in the case of stronger 
intersite coupling (figEF)- 




0.5 

(J> Ik (rad) 



FIG. 6: Effect of the strength of the intersite coupling J^d on the 
microscopic properties (as in figEJ in equilibrium (a,b,c,d) and at 
V = IV 

The graph in figd describes the distribution of sta- 
bility patterns of the three MMs in the constriction as 
a function of the exchange parameters. It is based on 
a study of the DW migration barrier for different ex- 
change parameters. We can distinguish three different 
regimes depending on the values of the exchange parame- 
ters: (1) "magnetostatic regime" associated with the pre- 
sense of 2 stable magnetic states, for which fa = ±7r/2; 
(2) "mixed regime": 4 stable configurations, 2 in each 
half- plane for which < fa < ir/2 and tt/2 < fa < 
7r; (3) "current-driven regime" : 8 stable configurations, 
namely all the collinear MM alignments (0, 0, 0), (ir, ir, ir), 
(0,0, tt), (O.tt.tt), (tt,0,0), (tt,^,0), (0,tt,0), (ttAtt) 
ffig llO|) . The last case is confirmed by full torque mini- 
mizations at various initial conditions. It should be noted 
that stability of the two unipolarized collinear configura- 
tions is also found above the boundary on figured but 
they are then only accessible once the system is trapped 
within very narrow regions of {fa space. Our main ob- 
servation is that the collinear configurations are the only 
form of stability of the magnetic chain in the contact for a 
range of exchange parameters around the "realistic" val- 
ues, defined earlier. Thus, even though our calculations 
in the previous and in the next section use Jdd = 50 meV 
and J = 1 eV, the qualitative features of the results may 
be expected to hold for a range of values of the parame- 
ters, with Jdd in the region of tens of meV and J > 0.5 eV. 

It is observed that the bias, driving a spin-polarized 
current, is able to distort the boundaries between these 
three regions significantly (figEI dashed lines refer to 
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FIG. 7: Diagram showing the three regions, with 2, 4 and 8 stable 
magnetic configurations respectively (see text), in a J-Jdd cu t of 
the parameter space. The boundaries between the three regimes are 
calculated from the properties of the (0, 0, 7r) — + (0, n, n) transition 
(see inset). The dashed lines correspond to voltages of 1 and 2 V. 



V = 1,2 V). The effect of the bias on the migration 
energy-barrier profile has been studied for different val- 
ues of the exchange parameters and reported for two rep- 
resentative cases J = 1 eV and J = 2.5 eV (figEJ- The 
overall observation is that for J < 1.5 eV (the "realis- 
tic" regime) increasing bias (current) enforces the bar- 
rier, while for J > 2.5 eV the bias suppresses the barrier. 
For intermediate values of J the barrier shows a non- 
monotonic behavior with bias. 
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FIG. 8: The biased energy barriers for the transition (0, 0, 7r) — > 
(0, 7T, tt) at different voltages V=0,1,2V, J = l.OeV (left) and J = 
2.5 eV (right) and for three values of Jdd (one in each of the ranges, 
discussed in the text): (a) Jdd = 12meV; (b) Jdd = 0.25 eV; (c) 
J dd = 0.3 eV; (d) J dd = leV; (e) J dd = 2eV; (f) J dd = 3eV. 



D. Current- voltage characteristics 

The current-voltage characteristics of the system at all 
the different stable alignments of the MM in the chain 
are presented in figure El All I-V curves are symmetric 
about the origin, due to the time reversal symmetry. The 
eight possible stable magnetic states are 4x2 degenerate 
with respect to the current, which, as discussed above, is 
due to the specific spatial symmetry of the atomic con- 
tact. The slopes of the I-V curves cannot be directly re- 
lated to the number of DWs in the constriction (three for 
the {0,tt,0} and {tt,0,0}; one for {0,0,0} and {0,0, tt} 
states). Nevertheless, the least steep curve docs corre- 
spond to the highest number of domain walls present 
within the chain. As soon as the DW migrates from the 
chain towards the leads (as in the {0,7r,0} to {tt, 0,0} 
transition) the conductance increases. We conclude that 
microscopic magnetization reversals in the constriction 
could be causing massive current variations (of up to 
50%) at a given bias. At a given finite temperature 
this would result in a random telegraph noise in conduc- 
tance measurements and such effects have been observed 
exp er iment ally£i2ii2S. 

Figure ITUI represents the work for successive 1 MM ro- 
tations. The sequence of transitions goes through every 
stable magnetic configuration once, and returns to the 
initial state. In equilibrium the depths of the wells in 
this graph correspond to the relative energies of our sys- 
tem in various stable magnetic states with respect to the 
initial one. Thus (0,7r,0) and (tt, 0, tt), which have two 
abrupt DW within the chain and are the least conduct- 
ing states, are found to be the most stable among the 
collinear alignments. It is observed that the external bias 
has a non-trivial effect on the effective energy barriers 
for these transitions. The total work for the closed loop 
cancels out for any bias. This, however, is not an indi- 
cation of conservativeness of the current-induced torques 
but is rather an artefact of the specific symmetries 
in this certain closed path, which includes all the states 
and can be decomposed into two branches going through 
identical states in opposite direction. 

There are 12 = (5 x 2 + 2) one-MM-rotation tran- 
sitions between different pairs of stable collinear align- 
ments, out of which we distinguish 5 + 2 = 7 different 
transitions. The average activation barrier at equilibrium 
is 71 meV with a variance of 36 meV and it depends non- 
monotonically on the bias: 65.6 meV at IV and 68meV 
at 2 V. These values for the activation barriers suggest 
switching frequencies, and hence random telegraph noise 
in the current, in the microwave range at room temper- 
ature. 



E. Are the torques conservative? 

The question if, and under what conditions, forces 
under steady-state current are conservative remains 
an open fundamental problem in the theory of 
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FIG. 9: The current- voltage characteristics at different stable 
collincar alignments for the MM in the chain. The inset repre- 
sents the correspondent spin-polarization of the net current kj = 
(7f ~ I [)/(!■{ calculated for each state. It is a non-monotonic 

function of the bias. 
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FIG. 10: Work for successive transitions between the 8 stable 
magnetic configurations at different bias V = 0, 1,2 V. 



transport^. A thermodynamic formulation of forces un- 
der non-equilibrium steady-state conditions, proposed in 
reference^, leads to the explicit identification of a ther- 
modynamic potential for electro-migrationS*. However, 
as a consequence of the infinite nature of open-boundary 
systems, this potential involves a conditionally conver- 
gent real-space summation. If the sequence of terms in 
this summation remains invariant along a given path in 
the configuration space of the system, then along that 
path, current-induced generalized forces are rigorously 
expressible as gradients of a scalar potential and are 
therefore conservative. The possibility remains open, 
however, that the order of terms in the conditionally con- 
vergent sum may change, as specific points, or manifolds, 
in configuration space are traversed^. This constitutes 
an effective breakdown of the Born-Oppenheimer approx- 
imation, with the consequence that paths that span such 
points are non-conservative2a. 

We now carry out a numerical test to see whether 
or not the torques in equation l|16[) are conservative in 
the present current-carrying system. The work for a set 
of transitions between collinear MM configurations, per- 
formed by rotation of a single MM, which form a closed- 
loop, is calculated for different voltages. The full work 



for three different loops of four consecutive transitions as 
a function of the applied voltage is presented in figure ITT1 
A significant variation of the closed-loop work with bias 
is observed. 

In order to resolve the numerical error we have per- 
formed a series of tests with different levels of accuracy. 
We recognize several sources of numerical error: (1) The 
level of convergence of the density matrix Sp. (2) The 
fineness of the energy mesh for the contour integration 
in the complex plane SE. (3) The level of torque relax- 
ation ST. (4) The angular mesh for the torque integra- 
tion which results into the work. As the torque in the 
current-driven regime is a very smooth function of the 
reaction coordinate (figHt), we have found (4) insignifi- 
cant for the value of the integral Q. The effect of the 
rest of the accuracy parameters on the closed loop work 
(see the loop {iririr — ttttO — 07r0 — Qirir — iririr} in fig lll|) 
is summarized in tabled 



(Sp,SE,ST), % 


Wov 


W0.5V WlV 


Wi. 5V 


W 2 v 


(100,100,100) 
(100,20,100) 
(1,100,2) 


0.0001 
0.0001 
0.0406 


-1.407 -7.454 
-1.290 -7.361 
-1.300 -7.356 


0.601 
-0.036 
0.662 


19.81 
19.79 
19.88 


TABLE I: The work (in mc 


V) for the loop {tttttt 


— 7T7rO — 


OttO - 



Otttt — tttttt} as function of the accuracy parameters (in relative 
units). The value of the bias V is given as a subscript. 



The results in table[I]suggest that the effect of nonzero 
work for a closed-loop set of transitions is not substan- 
tially affected by variations of 1-2 orders of magnitude 
about the chosen level of accuracy. Thus we see that 
in the present case, along the selected closed paths, we 
have an explicit example of non-conservative generalized 
non-equilibrium forces. 
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FIG. 11: Dependence of closed-loop work on the voltage 
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IV. CONCLUSIONS 

In conclusion, wc have proposed a microscopic 
quantum-classical approach for computing the current- 
induced torques on the local magnetization in ferromag- 
netic point contacts under bias. Our method employs an 
s-d model for the electronic structure and NEGF tech- 
nique for describing the electronic transport. The direc- 
tion of the local MMs are mapped onto classical degrees 
of freedom. Wc apply this method to a specific atomic 
structure, which consists of a monoatomic chain, bridg- 
ing over two semi-infinite leads with opposite magnetiza- 
tions, so that at least one magnetic DW is formed within 
the constriction. We investigate the stability of various 
magnetic configurations, involving multiple DWs, and 
the effect of bias driving a spin-polarized current, on the 
energy-barrier for the DW migration. For realistic values 
of exchange parameters only the collinear MM arrange- 
ments are stable. These configurations carry different 
(by up to 50 %) net currents, and the average activa- 
tion barrier for transitions is about 65-70 meV with vari- 
ance of 20-40 meV, depending on bias. Therefore ran- 
dom telegraph noise in current with significant amplitude 
could be related to thermally-activated MM rearrange- 
ments within the constriction. We have also found that 
geometrical asymmetries in the atomic structures (which 
could be induced by the current^) affect the symmetry of 
the activation barrier for DW migration, pinning it to a 



preferential spatial position, in which case the structure 
becomes less conducting. 

The observation that the collinear MM alignments 
are the only stable magnetic states and the fact that 
direct intersite interaction is not able to induce flips 
in the neighboring MMs as one MM in the chain is 
quasi-statically rotated, has enabled us to calculate the 
work for series of successive magnetic rearrangements of 
MMs, involving single-MM rotations. Thus, we could 
address numerically the long-standing question of the 
conservativeness of the current-induced forces (torques, 
in our case) in open-boundary non-equilibrium system. 
We have found numerical evidence that the work for 
various closed-loop paths is not zero, but varies non- 
monotonically as the system is driven away from equi- 
librium. Hence generalized current-induced forces in our 
present system are not conservative, at least in the sec- 
tion of the configuration space spanned by the present 
calculations. 
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